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Abstract 

The influence of electron-electron (e-e) interactions on the transmission 
through a quantum dot is investigated numerically for the Coulomb blockade 
regime. For vanishing magnetic fields, the conductance peak height statis- 
tics is found to be independent of the interactions strength. It is identical to 
the statistics predicted by constant interaction single electron random matrix 
theory and agrees well with recent experiments. However, in contrast to these 
random matrix theories, our calculations reproduces the reduced sensitivity to 
magnetic flux observed in many experiments. The relevant physics is traced to 
the short range Coulomb correlations providing thus a unified explanation for 
the transmission statistics as well as for the large conductance peak spacing 

fluctuations observed in other experiments. 
PACS numbers: 71.55.Jv,73.20.Dx,71.27.+a 
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Since the discovery of the Coulomb blockade phenomenon, most tunneling experiments 
through a quantum dot were interpreted within the constant interaction model (CI) 
In that approximation, the ground state energy of a quantum dot populated by N electrons 
is expressed as En = + X)£Li Vi where C is the dots constant (or slowly varying) 
capacitance, and rji are the single particle energies. Evidently, in that model only the 
long range Coulomb interaction is taken into account while the short range correlations 
are neglected. The fine ground state properties are hence determined by the single particle 
states which for a disordered or chaotic dot display a random matrix theory (RMT) statistics. 
Although the CI model is very appealing in its simplicity, it was recently proved wrong in 
predicting the the distribution of the conductance peak spacings [§-0], as well as the large 
peak spacing fluctuations found in some of the experiment . While the CI model predicts 
a RMT type of ground state statistics with a characteristic energy scale, A (average single 
particle level spacing), the experiments find a different type of statistics and (at least for 
the experiments in Ref. 0,0) considerably larger fluctuations which seem to be independent 
of A. Moreover, the ground state energy turns out to be relatively insensitive to magnetic 
flux ||] and application of one quantum flux unit, (f> = hc/e, through the dot hardly affect 
it [0,0. This insensitivity is again in contrast with the CI model since the single particle 
states, and hence the ground state energy in that model, are expected to fluctuate on a flux 
scale smaller than one quantum flux unit HH. 

A point of view similar to the CI one was also taken for the calculation of the Coulomb 



blockade conductance peak height statistics ||iO|-[T3|| . Since RMT was assumed for the single 



particle states, a modified Porter-Thomas distribution was obtained for the dimensionless 
transmission a = 2TlTr/T(T) 

r 2 



P s =o(«) = \ — e 
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PbM<*) = 4a[K {2a) + K^a)}^ 2 *. (1) 

Here, ri(T^) are the tunneling rates from the left (right) lead, r = Tl + Tr, (. . .) denotes 
an average over different peaks or disorder realizations and K ,Ki are the modified Bessel 



functions. 

Subsequent experiments |||7| reported partial agreement with these calculations and one 
was therefore facing the following dilemma: while the conductance peak spacing fluctuations 
can be explained by short range Coulomb correlations || (see also recent papers by Koulakov 
et al. jnj and Blanter et al. JI5|), the conductance peak height statistics roughly agrees with 



a model that totally neglects these correlations. It is hard to reconcile such two different 
pieces of physics for two facets of the same phenomenon and in the present manuscript 
we show indeed that the observed peak height statistics may result from fluctuations in 
the short range Coulomb correlations rather than the single electron RMT physics (i.e., 



wave functions with no correlations) utilized in refs. [|T0| p~3|] . Strong support in favor of the 
Coulomb correlations type of physics comes from the insensitivity to magnetic flux observed 
in both experiments |||7]]. The auto correlation function between the height of a given peak 
at two different values of the magnetic flux, 0, 

(fo(£)fo( £i A£ L; fo = 

is found experimentally to decay on flux scales larger than predicted by CI single electron 
RMT. This reduced sensitivity to flux is clearly manifested in the numerical calculation pre- 
sented below that take the interactions into account. Curiously, the inclusion of interactions 
does not change the peak hight distribution predicted by single electron RMT and observed 
by the experiments. We are therefore able to propose a unified explanation for both facets of 
the Coulomb blockade phenomenon, namely, the conductance peak spacing fluctuations and 
the peak height distribution. They both may originate from fluctuations in the Coulomb 
interaction rather than single particle physics 

The height g ma x of a conductance peak is given by 

e 2 / 7r 



The tunneling rates may be formulated in a tight-binding many particle language as 

r L (R) = \t L{ R)\ 2 \ E (^N+MJ^n)^, (4) 

k,j<E[L(R)] 



where turn is the barrier transmission which is assumed to depend only weakly on energy, 
\l/jv is the N particle ground state wave function in the dot, a k i is the fermionic creation 
operator at site (k,j), and summation is performed on sites k,j 6 [L(R)] i.e, sites adjacent 
to the left (right) lead. This is a straight forward adaptation of the definition of r^(^) given 
in Ref . |Tl| . For the non- interacting case, assuming statistically identical independent single 



channel leads (for example, point contacts |TT]]), RMT predicts the distribution depicted in 
Eq. © 0]. 

We calculate the tunneling rates for a system of interacting electrons modeled by 

a tight-binding Hamiltonian. We choose a 2D cylindrical geometry of circumference L x and 
height L y . This particular geometry is very convenient for the study of the influence of a 
magnetic flux threading the cylinder in the y direction. A radial magnetic field could 
also be applied, but for a field equivalent to one quantum flux unit through the system, 
Landau bands appear. Since this is not the situation in the experiment we prefer to apply 
a threading flux only. The Hamiltonian is given by: 

H = J2 e *J a h a *J - v Y.^^ i27[ ^l^)s/L x )al j+l a Kj + a\ +l j a Kj + h.c) + H int , (5) 

k,j k,j 

where e k j is the energy of a site (k,j), chosen randomly between — W/2 and W/2 with 
uniform probability, V is a. constant hopping matrix element, and s is the lattice unit. The 
interaction Hamiltonian is given by: 

H int = u y, ^ ak fh" l ; p (6) 

k,j>l,p \ r k,j r l,p\/ S 

where U = e 2 /s. The distance \r k j — fi tP \/ s = (min{(fc — I) 2 , (L x /s — (k — Z)) 2 } + (j — p) 2 ) 1 / 2 . 
The interaction term represents Coulomb interaction between electrons confined to a 2D 
cylinder embedded in a 3D space. 

We consider a 4 x 6 cylinder with M = 24 sites and iV = 3 or iV = 4 electrons. The 
size of the many body Hilbert space is m = The m x m Hamiltonian matrix is 

numerically diagonalized and the ground-state eigenvectors are obtained. The strength 
of e-e interactions, U, is varied between — 22V. The disorder strength is set to W = 3V 



in order to assure RMT behavior for the non-interacting case. For each value of U, the 
results are averaged over 500 different realizations of disorder. The left lead is attached to 
the (1, 1) site and the right one to the (4, 6) site (point contacts). Assuming a fixed barrier 
transmission |ti(_R)| 2 = 1, Tl(R) is calculated using Eq. (f|). 

The average tunneling rate, its root mean square value, and the normalized fluctuations 
are presented in Fig |l] as function of the interaction strength U for B = and B ^ 0. Using 
the expression of ^l(r) i n terms of Green functions given by Zyuzin and Spivak Hlq 



one can 



use a diagrammatic summation similar to the one used in Fig. 1(a) of Ref. |18| for dN/dfi 
to obtain the RPA predictions for (Tl(r) (U)). For small values of U the RPA prediction [I 



(Fl(r)(0))/(T L (r)(U)) ~ 1 + (kL/tt) ~ 1 + 0.34U (where k = S d ue 2 , S d = 2vr for an infinite 
system and 2.50 for a 4 x 6 lattice [18|, v = (AL 2 ) -1 and L 2 = L x L y ) is followed for both 
B = and B ^ 0. For larger values of U, the tunneling rate, T L ^(U) is reduced below 
the RPA value. It is expected that the short range order induced by the interactions indeed 
reduce the overlap between a^\^^) and l^jv+i}- 

The fluctuations in the tunneling rate depicted in Fig. |I|b, also decrease as function 
of the interaction strength. In the RPA regime the fluctuations may be calculated by a 
diagrammatic expansion similar to the one used in Fig. 1(b) of Ref. [[TJ| for 5 2 dN/dfi. 



As in Refs. fllqjlql a cutoff must be used in order to avoid divergences in the non- 
interacting case, but the influence of the interaction is cutoff independent resulting in 
((5 2 r L(i?) (0))/(5 2 r L(i?) (f/))) 1 /2 = (l + ^L/vr)) 2 , while for stronger interactions (5 2 T L{K) (U)) 
is strongly suppressed. 

To make connection to real samples we use the ratio between the average inter-particle 
Coulomb interaction and the Fermi energy r s = 1/ y/nna b (where n is the electronic density 



and Og is the Bohr radius) corresponding to r s ~ yn/Q(U/2V) for iV = 4, M = 24. For 
all the above quantities there is a clear borderline around U = 2 — AV or r s ~ 1. At low 
r s values (high densities) the tunneling rate agrees well with RPA calculations, while for 
stronger interactions the results are qualitatively different. Identical behavior was observed 
for conductance peak spacings fluctuations discussed in Ref. @. In both cases the appearance 



of short range correlations at r s > 1 lead to a failure of RPA. It is important to bear in 
mind that for real systems the density in the leads n = 2 — 3.5 x 10 11 cm -2, while in the 
dot the density is probably lower, thus in all experimental systems r s ~ 1 — 2. They hence 
correspond to a regime where RPA no longer holds. 

The full probability distribution of the dimensionless parameter a for different values 
of the interaction is shown in Fig. ||]. For B = the distributions for all values of U are 
reasonably close to the RMT prediction, Eq. (fj). Thus, moderate changes in the second 
moment of r^(^) (Fig. |l]c) hardly influence the distribution of a. This is in good agreement 
with the experimental data |||7|. Also in the presence of a magnetic field, the interaction 
strength has no major effect. Nevertheless, it is interesting to note that the dip at small 
values of a predicted by RMT (Eq. (P) seen for small values of U (actually the dip is even 
larger than predicted which is an artifact of the small system size) disappears for higher 
values of U . The disappearance of the dip at higher interaction values is the result of the 
reduced sensitivity of the system to the magnetic field. 

The auto-correlation, Eq. (|2[), between the height of the i-th peak at different values of 
magnetic fields is given according to RMT P~3,[I3| for the GUE ensemble by 



C(<j>,A<j>) 



(7) 



Here <ft = BA is the magnetic flux through a dot of area A, and <ft c = (po/^gJC, where for 
a diffusive dot g is the dimensionless conductance g = E c / A (E c is the Thouless energy) 
and K, = 1. For a ballistic dot, E c = vp/y/A, and the geometrical factor, for the case of a 
flux line threading the dot K is of order of unity [|13|]. This corresponds to <p c ~ 0.1 O f° r 
the experimental setup in Ref. 0, while the experimental value is C ~ 0o Hl9 ]- Recently 



Alhassid ]2(J pointed out that for a uniform magnetic field in the dot /C is much smaller 
corresponding to C ~ O.50o, which is still significantly lower than the experimental value. 

The numerically calculated C(<j), A<j)) for cf) = O.20o for different values of U is presented 
in Fig |3|. Eq. (|7|) describes the U = behavior for small values of A0 quite well. From the 
value of 0c = O.50o one obtains g = 4 which is reasonable. 
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Since the <fi dependent part in the diagrammatic calculation (which is similar to the 



calculation of the fluctuations p!6|JT8|) does not depend on the interaction, one expects to 



first order no changes in Eq. ([!]) from RMT calculations [|T/]]. Indeed, for small U there is 
only a weak influence of the interaction strength on the auto-correlation function. At larger 
values of A<f>, C ((f), A0) seems to be more sensitive to U and in the regime V < U < AV it 
is hard to fit it by a particular functional form. For U > AV (r s > y2) the auto-correlation 
function can be fitted again to the functional form of Eq. (0) but the correlation flux is 
enhanced. For example, for U = 10V <p c = 1.750o- Thus, although strong interactions reduce 
the average conductance peak height, they hardly influence the conductance distribution, 
again in agreement with both experiments ||[7|]. The main effect is a reduced sensitivity 
to flux which is also manifested in the weak dependence of the conductance peak spacing 
fluctuations on magnetic field. It is worth mentioning that tunneling experiments through 
excited states of heavily doped GaAs dots 0, where r s < 1, find the RMT correlation flux. 
The numerical calculation gives the two particle correlation function as well ||18|| . It turns out 
that the point where the conductance starts to deviate from RPA (r s ~ 1) is accompanied 
by the appearance of short range correlations, the same correlations that led to the large 
conductance peak spacing fluctuations are here responsible for the transmission statistics. 
We emphasize that Wigner crystallization occurs at much stronger interactions. 

In conclusion, the appearance of short range electron correlations which are the result 
of the low electronic densities in the measured quantum dots explains not only the large 
conductance peak spacing fluctuations, but also other features pertaining to their transport 
properties. First, it has been numerically shown that although the conductance peak hight 
distribution is well described by the single electron RMT, this distribution is valid also at 
the low density (strong interaction) regime. The strong interaction model also agrees with 
the results of the experiment in the presence of a magnetic field, which stands at odds with 
the predictions of the single electron RMT. Most importantly, the auto-correlation between 
peak heights exhibits an enhanced characteristic magnetic field for r s > 1 which is in good 
agreement with the experimental observations Pj7],[n|]. Thus, the fact that the experiments 



are performed at low densities (r s > 1) for which the RPA theory no longer holds explains 
many of the puzzling features exhibited by these systems. A unified explanation for the 
large conductance peak spacing fluctuations and the peak height characteristics is hence 
provided. A better understanding of temperature effects, and a full analytic treatment of 
the short range correlations is still needed. 
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FIG. 1. The influence of the e-e interaction strength on the tunneling rate ^l(r)- Circles 
correspond to <j> = 0, while squares to <p = O.40o- The lines represent the prediction of an RPA 
theory, (a) The average tunneling rate, (b) the variance. 
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FIG. 2. The probability distribution of the dimensionless parameter a corresponding to the 
conductance peak height g max for different values of the e-e interaction strength, (a) 4> = 0, (b) 
4> = O.40Q. The lines correspond to the RMT predictions for B = (GOE) and B / (GUE). 
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FIG. 3. The auto-correlation function for different values of the e-e interaction strength. The 
full line corresponds to Eq. ([?]) with (j) c = 0.5(/>o, while the dashed line to Eq. (|7|) with <f> c = 1.75(/>o. 
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